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Abstract: We present a family of one- dimensional cellular automata mod- 
eling the diffusion of an innovation in a population. Starting from simple 
deterministic rules, we construct models parameterized by the interaction 
range and exhibiting a second-order phase transition. We show that the 
number of individuals who eventually keep adopting the innovation strongly 
depends on connectivity between individuals. 



1 Introduction 

Diffusion phenomena in social systems such as spread of news, rumors or in- 
novations have been extensively studied for the past three decades by social 
scientists, geographers, economists, as well as management and marketing 
scholars. Traditionally, ordinary differential equations have been used to 
model these phenomena, beginning with the Bass model (|Bass, 1969|) and 



ending with sophisticated models that take into account learning, risk aver- 
sion, nature of innovation, etc. ([Mahaian et al., 1990|; Mahajan and Peterson 



1985|; Rogers, 1995| ). Models incorporating space and spatial distribution of 



individuals have been also proposed, although most research in this field has 
been directed to refining of the discrete Hagerstrand models (jHagerstrand 



1952, 1965|) and to constructing partial differential equations similar to dif- 



fusion equations known to physicists ( |Haynes et al., 1977] ) 



Diffusion of innovations (we will use this term in a general sense, mean- 
ing also news, rumors, new products, etc.) is usually defined as a process by 
which the innovation "is communicated through certain channels over time 
among the members of a social system" ( [Rogers, 1995 ). In most cases, these 
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communication channels have a rather short range, i.e. in our decisions we 
are heavily influenced by our friends, family, coworkers, but not that much 
by unknown people in distant cities. This local nature of social interactions 
makes cellular automata (CA) a well-adapted tool in modeling diffusion phe- 
nomena. In fact, epidemic models formulated in terms of automata networks 
have been successfully constructed in recent years ( [Boccara and Cheong 
19921, 19931; Boccara et al, 1994|) . 



2 Simple deterministic models 

We will construct models of diffusion of innovations based on elementary 
(radius-1) cellular automata ([Wolfram, 1986|) . If s{i,t) denotes the state of 



lattice site i at time t, the function {i, t) i— >■ s(z, t) is a mapping from Z x N to 
Z2 = {0, 1}. Given a function / : {0, Ij^'^+i 1-^ {0, 1}, the discrete dynamical 
system 

s{i, t + l) = f{s{i - r, t), - r + 1, t), . . . , s(z + r, t)) (1) 

is called a cellular automaton (CA) of radius r, and / is called its local 
(function) rule. 

In the simplest version of our model, the sites of an inflnite lattice are all 
occupied by individuals. The individuals are of two different types: adopters 
(^s{i,t) = ij and neutrals (^s{i,t) = 0^. Moreover, we wil assume that an 
individual can get information only from it two nearest neighbors. To simplify 
our model, we will also assume that once an individual becomes an adopter, 
he remains an adopter, i.e. his state cannot change. This condition flxes 
four entries in the rule table: {0, 1, 0} 1, {0, 1, 1} ^ 1, {1, 1, 0} — 1 and 
{1, 1, 1} 1, and since the information comes from nearest neighbors, if 
both are neutral, then the individual will stay neutral, i.e., {0,0,0} 0. 
This leaves three entries in the rule table to be determined, and this can be 
done in 2^ = 8 ways, as shown in Table |I|. The flrst Rule 204 (for rule codes 
cf. Wolfram 86) listed in this table is just the identity, and it will be excluded 
from further considerations. Rules 220 and 252 can be obtained respectively 
from 206 and 238 by spatial reflection, therefore they will be excluded too. 
This leaves us with flve distinct rules. 

Rule 254: An individual adopts if, at least, one of his neighbors is an 
adopter. 
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Table 1: Eight possible elementary rules for diffusion of an innovation. 

Rule 238: An individual adopts only if his right neighbor is an adopter. 

Rule 222: An individual adopts only if exactly one of his neighbors is an 
adopter. 

Rule 206: An individual adopts only if his right neighbor is an adopter and 
its left neighbor is neutral. 

Rule 236: An individual adopts only if both his neighbors are adopters. 

We will now demonstrate that in all five cases, the density of adopters p{t) 
at time t can be exactly computed, assuming that we start from a disordered 
initial configuration with p(0) = po- 

(i) In order to understand the dynamics of Rule 254, we can view the 
initial configuration as clusters of ones separated by clusters of zeros. Only 
neutral sites adjacent to a cluster of ones change their state, while all other 
neutral sites remain neutral, as shown in the example below (sites that will 
change their state to 1 in the next time step are underlined): 
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This implies that the length / of every cluster of zeros decreases by two every 
time step, i.e. 

M{l,t + 1)^ M{l + 2,t), (2) 

where M{1, t) denotes a number of clusters of zeros of size s per site at time 
t, and therefore 

M{l,t) ^ M{l + 2t,Q). (3) 
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For a random initial configuration with initial density the cluster density 
is given by M{1, 0) = (1 — Po) Vo> hence 



\l+2t 2 

Po- 



M(/,t) = (l-po)' 
The density of zeros at time t, denoted by r}{t), is given by 



77(t) = ^/M(/,t) = (l-po) 
1=1 



2t+l 



and finally the density of ones p{t) is equal to 

pit)^l-rjit)^l-il-por^\ 



(4) 



(5) 



(6) 



(ii) For Rule 238, the derivation is similar, except that now every cluster 
of zeros decreases by one every time step, which means that we have to 
replace 2t by t in the previous result, i.e. 



p(t) = l-(l-po) 



t+i 



(7) 



(iii) In Rule 222, individuals "dislike overcrowding", and they do not 
become adopters if their two neighbors are already adopters, adopt if both 
neighbors adopted, i.e. {1, 0, 1} — > 0. As a consequence, clusters of size / > 1 
decrease their size by two units every time step, while clusters consisting of 
a single isolated zero (/ = 1) do not change their size. Clusters of size 1 are 
thus created from clusters of size 3, as well as those of size 1. For the density 
of ones, this yields (using a similar reasoning as before) 



p{t) = l-pl{l-po)-po 



2-po 



2 , Po(l - Pof 

Po ' 



Po 



[i-por'^\ (8) 



and 



lim p{t) = p{oo) = 1 - pI{1 - Po) - Po 

t— »oo 



:i - PoY 



2 - Po 

(iv) Rule 206 is just like Rule 222, except that clusters of zeros decrease 
their length by one. Therefore, we obtain 



and p(oo) 



p(t) = l-po(l-po)-(l-po)*+', 
1 -po(l -po). 



(10) 
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(v) Rule 236 is the simplest because, after the first iteration, all clusters of 
zeros of size 1 disappear while all other clusters remain unchanged, therefore 

p(t) = l-(l-po)'(l + Po). (11) 

In summary, in all cases (except for Rule 236), the density of ones ap- 
proaches the fixed point p{oo) exponentially. In a real social system, however, 
this is not the case. The density of adopters usually follows an S-shaped or 
logistic curve. The model discussed in the next section eliminates this short- 
coming. 



3 Probabilistic model 

In order to generalize the simple model discussed previously, consider a 2- 
state probabilistic cellular automaton, with a dynamics such that s{i,t + 1) 
depends on s{i,t) and (j{i,t), where 

oo 

ait,t)= J2 sit + n,t)p{n), (12) 

n=— oo 

and p is a nonnegative function satisfying 

oo 

E P(^) = 1- (13) 

n=— oo 

(deterministic automata networks of this type have been studied in details 
by Poccara et al., 19971) . 

Our generalized model is defined as follows. At every time step, a neutral 
individual located at site i at time t can become an adopter at time t + 1 with 
a probability depending on a{i, t) (in this section we will simply assume that 
this probability is equal to a{i,t)). As before, once an individual becomes 
an adopter, he remains an adopter, i.e. his state cannot change. 

The model can be viewed as a probabilistic cellular automaton with the 
probability distribution 

P(s(z,t + 1)=0) = {l-s{i,t)){l-a{i,t)) (14) 
P(s(i,t + 1) = l) = I - (l - s{i,t)){l - a{i,t)) (15) 
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The transition probability Pb^a is defined as 

Pb^a = P[s{i,t+l) = h\s{i,t) = a), (16) 

and represents the probabihty for a given site of changing its state from a 
to h in one time step. In our case, the transition probabihty matrix has the 
form 



Po^o Po^i ^ l-a{i,t) 
Pi^o Pi^i \ [ (r{i,t) 1 

As a first approximation we we consider (T{i,t) defined by 



-1 



R 



f^(^,^) = ^( J2 s{i + n,t) + J2sii + n,t) 



(17) 



;i8) 



(j{i,t) is then the local density of adopters at time t over the 2R closest 
neighbors of site i. This choice of a, although somewhat simplistic, captures 
some essential features of a real social system: the number of infiuential 
neighbors is finite and these neighbors are all located within a certain finite 
radius R. Opinions of all neighbors have equal weight here, which is maybe 
not realistic, but good enough as a first approximation. Let p(t) be the global 
density of adopters at time t (i.e. number of adopters per lattice site), and 
ri{t) = 1 — p{t). Since Pi^i = 1, the number of adopters increases with time, 
and limj^oop(^) = 1- If we start with a small initial density of randomly 
distributed adopters po, pit) follows a characteristic S-shaped curve, typical 
of many growth processes. The curve becomes steeper when R increases, and 
if R is large enough it takes only a few time steps to reach a high density 
of ones (e.g. 0.99). Figure |l| shows some examples of curves obtained in a 
computer experiment with a lattice size equal to 10^ and po = 0.02. 



4 Exact Results 

Average densities p{t+l) and r]{t + l) can be obtained from previous densities 
pit) and rjit) using the transition probability matrix ( () denotes a spatial 
average) 

" r/(t+l) 
_p{t + l) 

hence 

Pit + 1) = p(t)(Pi^i) + (1 - p(t))(Pi^o), (20) 



(Po^o) (Po^i) 
(Pi^o) (Pi^i) 



vit) 

Pit) 



(19) 
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Figure 1: Density of adopters (vertical axis) as a function of time (horizontal 
axis) for several values of radius R. 
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and using (|T^ we have 



p(t + l) = pit) + {l-pit)){ait,t)). (21) 

This difference equation can be solved in two special cases, R = 1 and R = oo. 

If -R = 1, only three possible values of local density a are allowed: 0, | and 
1. The initial configuration can be viewed as clusters of ones separated by 
clusters of zeros. Only neutral sites adjacent to a cluster of ones can change 
their state, and they will become adopters with probability |. All other 
neutral sites have a local density equal to zero, therefore they will remain 
neutral. This implies that the length of a cluster of zeros will on the average 
decrease by one every time step, just like in the case of Rule 238. The density 
of adopters, therefore, is, as for Rule 238, given by 

Pit) = 1 - (1 - poY^'. (22) 

Using the above expression, we obtain 

p(t + l) = p(t) + (l-p(t))po. (23) 

Comparing with ( pT]) this yields {a{i,t)) = po, i.e. the average probability 
that a neutral individual adopts the innovation is time-independent and equal 
to the initial density of adopters, which was a priori not obvious. 

When i? — > oo, the local density of ones becomes equivalent to the global 
density, thus in (PT|) we can replace {a{i,t)) by p(t). Hence 

p(t + l) = p(t) + (l-p(t))p(t), (24) 

or 

l-p(^ + l) = (l-p(t))^ (25) 

that is 

p{t) = l-il-pof. (26) 

Note that this case corresponds to the mean-field approximation, as we ne- 
glect all spatial correlations and replace the local density by the global one. 

For 1 < < oo, we may assume that the density of adopters can be 
written in the form 

p(t) = l-(l-po)^(*'^), (27) 
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where / is a certain function satisfying f(t, 1) = t + 1 and f{t, oo) = 2*. 
Computer simulations suggest that for a finite R, f(t, R) becomes asymptot- 
ically linear (when t — > oo), and the slope of the asymptote increases with 
R. Moreover, for large R, f satisfies the following approximate equation 

/(t + fc,2^i?) = 2V(t,/2), (28) 

where k is a positive integer. Detailed discussion of f(t, R) and its properties 
can be found in [buks and Boccara (1996| ). 



5 Generalization 



The model presented in the previous chapter was rather crude. One of its 
assumptions is that once the individual accepts the innovation, he will never 
change his mind. In practice, every technology or product has a finite life 
span. For some products, as TV sets, this time is relatively long, while for 
other items, like computer software, it is much shorter. One way to incor- 
porate this phenomenon in our model is to assume that at every time step, 
any adopter can drop the innovation with a given probability p. Therefore, 
the average time during which an individual is an adopter is 1/p (geometric 
distribution) . 

To make the model even more realistic, we will also assume that the 
adoption probability is not equal but proportional to the local density of 
adopters 

Pi^o = g^(^,t), (29) 
where q G [0, 1]. Hence, the new transition probability matrix is 



Prr 



Po^i 

Pl^l 



1 - qa{i,t) 
qa{i,t) 



P 

1 — p 



The difference equation for the average density of adopters 
p{t + 1) = p{t) (Pi^i) + (l - Pit)) (Pi^o) 



now becomes 



Pit + !) = (!- p)pit) + g(l - Pit)) (a(z, t)). 



(30) 



(31) 



(32) 



The previous model is recovered for p = and q = 1. When q = 0, the solu- 
tion is pit) = po(l —p)^- When R = 1, the model can also be considered as a 
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discrete version of the contact process, an irreversible lattice model involving 
nearest-neighbor interactions, used to study catalytic reactions (Parris, 1974 



Liggett, 1985| ). In the contact process, a particle desorbs spontaneously with 
rate p and adsorbs at a given unoccupied site at a rate proportional to the 
number of neighboring occupied sites. 

Although it is not possible to solve exactly equation ([3^), the general 
nature of the solution can be understood using the mean-field approximation 
(MFA), in which it is assumed that the average local density {a{i,t)) is the 
same as the global density p(t) (this is actually true for R = oo). Our 
equation becomes 

p(t + 1) = (1 - p)p{t) + qpit) (l - pit)) . (33) 

The substitution 



yields 

Xt+l 



Pit) = ^— ^a;(t) (34) 



p + q)x{t)(l-x{t)). (35) 



For R = oo, therefore, the dynamics of the model can be understood as an 
iteration of the logistic map Qx : x Xx{l — x) with X = 1 — p + q. Note 
that < A < 2, which excludes stable periodic points and chaos. 

Qx has always two fixed points a;'^^)(oo) = and a;*^^^(oo) = 1 — 1/A. Only 
one, however, is stable, depending on A, namely a;*^^''(oo) when A < 1 and 
x^'^\oo) when A > 1. In terms of p{t) we obtain 

/il^'^(^) I l—p/q otherwise. ^^^^ 

Since |(5^(xoo)| < 1 if A 7^ 1, the stable fixed point is hyperbolic (strongly 
attracting). At A = 1, however, it becomes nonhyperbolic (weakly attract- 
ing), and Qx exhibits a transcritical bifurcation (|Bardos and Bessis, 19801) 
with exchange of stability. 



6 Phase Transition 

The mean-field approximation discussed in the previous section becomes cor- 
rect only when R ^ 00. For small values of R, like in the basic model, 
strong correlations are created and substantial deviations from MFA can be 
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Figure 2: [p, q) phase diagram obtained using Monte Carlo simulations. Lines 
separating regions where p(oo) = and p(oo) ^ are shown for i? = 1 (o), 
R = 4 (□), and R = 16 (A). The dotted line represents the mean-field 
approximation. Lattice size is equal to 10000. 

expected. Figure H represents the (p, q) phase diagram for different values of 
R obtained in computer simulations. The smaller R is, the larger the devi- 
ation from MFA (dotted line) becomes. As we can see, the line separating 
p(oo) 7^ and p(oo) = shifts to the left (toward larger p) when R increases. 
This means that the connectivity between sites increases the robustness of 
the innovation, i.e. even if p is large, it can be compensated by a large R. 
To see it, consider a point on the (p, q) phase diagram which is located be- 
tween the R = 1 and the mean-field lines (e.g. (0.2,0.3)). If we plot p(oo) 
as a function of R, we obtain a bifurcation (or phase-transition) diagram, 
as shown in Figure For R < 6 the asymptotic density of adopters goest 
to zero, but if i? > 6, p{oo) becomes positive. Three distinct regions of the 
(p, q) phase space, therefore, can be distinguished: 

1. In the region bounded by the q axis and the R = 1 critical line, p{oo) > 
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Figure 3: Phase transition of p(oo) with R playing the role of the control 
parameter for p = 0.2 and q = 0.3. 

regardless of R. 

2. In the region bounded by the R = 1 and the mean-field critical lines, 
p(oo) is either equal to zero or positive, depending on the value of R. 
In general, p(oo) is positive for a sufficiently large R. 

3. In the region bounded by the mean-field critical line, p(oo) = regard- 
less of R. 



7 R = 1 critical line 

We studied the boundary between regions 1 and 2, i.e., the R = 1 critical 
line, using the local structure theory (LST) up to order 4 ( |Gutowitz, 1987 
[Gutowitz et al., 1987|) , as shown in Figure |[ 
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Figure 4: R = 1 critical line (□) and its local structure theory approxi- 
mations. Starting from the right, the consecutive solid lines correspond to 
n = 1, 2, 3 and 4. 



As we increase the order n of the approximation, we obtain an approx- 
imate critical line closer to the experimental one. This can be also seen if 
we plot p{oo) obtained from LST as a function p, keeping q constant. An 
example of such curves is shown in Figure ^, where six consecutive orders 
of LST are presented. For this particular value of p, the sixth-order critical 
g- value qc differs by about 11% from the numerical result, as shown below. 



n 


1 


2 


3 


4 


5 


6 


simulation 


(Ic 


0.200 


0.362 


0.430 


0.462 


0.479 


0.490 


0.549 



As mentioned earlier, for q = 1 and R = 1, our model can be consid- 
ered as a discrete realization of the contact process. For the contact pro- 
cess, theoretical studies as well as Monte Carlo simulations suggest that the 
phase transition occurs at Pc=0.3032 (Pickman and Jensen, 1991|; Jensen and 



Dickman, 1993|) . The table below compares this result with our numerical 
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Figure 5: Local structure approximation for p — 0.2. Solid lines represent 

consecutive orders of the local structure approximations, starting with n = 1 
(the leftmost curve), and ending with n = 6 (the rightmost solid curve). The 
dots represent results of computer simulations (p = 0.2, R = 1). 

simulations and local structure theory approximate results. 
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The discrepancy between discrete and continuous time versions of the contact 
process is understandable, since the synchronous dynamics of CA allows, 
for example, for one site to change its state from 1 to 0, and, at exactly 
the same time, for its neighbor to change its state from to 1, something 
that does not occur in the continuous time model. We should stress at this 
point that although there are many analytical techniques developed to treat 
stochastic processes such as the contact process, none of them, in general, 
can be easily translated to the language of cellular automata theory, and 
the synchronicity of updating is usually the main source of difficulties. For 
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example, the cluster approximation recently proposed for the contact process 
(|Ben-Naim and Krapivsky, 1996|) belong to this category. 

On the other hand, when all transition probabilities are small, a single 
site of the lattice is rarely changed, and during most of the iterations of the 
CA it will remain in the same state. The updating, therefore, becomes more 
and more similar to asynchronous updating, and as a consequence, we can 
expect that the model should behave almost like the time-continuous contact 
process. 



8 Conclusion 

We have studied a probabilistic cellular automata model for the spread of 
innovations. Our results emphasize the importance of the range of the inter- 
action between individuals. The innovation spreads faster when the range 
increases since increased connectivity between individuals reduces constrains 
on the exchange of information. Larger connectivity could be also achieved 
by increasing space dimensionality. 

The range of interaction R not only affects the growth rate, but in the 
region of parameter space {p, q) bounded by the R = 1 critical line and 
the mean-field line it can be a decisive factor for the asymptotic density of 
adopters p(cxo). In this region, if R is too small, p(cxo) = 0, while for a 
sufficiently large connectivity p{oo) > 0. 
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